\(\int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx\) [128]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [F]
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 20, antiderivative size = 85 \[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=-\frac {\cos ^2(a+b x)^{\frac {1-m}{2}} \csc ^2(a+b x) \operatorname {Hypergeometric2F1}\left (\frac {1-m}{2},\frac {1}{2} (-2+m),\frac {m}{2},\sin ^2(a+b x)\right ) \sec (a+b x) \sin ^m(2 a+2 b x)}{b (2-m)} \]

[Out]

-(cos(b*x+a)^2)^(1/2-1/2*m)*csc(b*x+a)^2*hypergeom([-1+1/2*m, 1/2-1/2*m],[1/2*m],sin(b*x+a)^2)*sec(b*x+a)*sin(
2*b*x+2*a)^m/b/(2-m)

Rubi [A] (verified)

Time = 0.09 (sec) , antiderivative size = 85, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.100, Rules used = {4395, 2657} \[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=-\frac {\csc ^2(a+b x) \sec (a+b x) \sin ^m(2 a+2 b x) \cos ^2(a+b x)^{\frac {1-m}{2}} \operatorname {Hypergeometric2F1}\left (\frac {1-m}{2},\frac {m-2}{2},\frac {m}{2},\sin ^2(a+b x)\right )}{b (2-m)} \]

[In]

Int[Csc[a + b*x]^3*Sin[2*a + 2*b*x]^m,x]

[Out]

-(((Cos[a + b*x]^2)^((1 - m)/2)*Csc[a + b*x]^2*Hypergeometric2F1[(1 - m)/2, (-2 + m)/2, m/2, Sin[a + b*x]^2]*S
ec[a + b*x]*Sin[2*a + 2*b*x]^m)/(b*(2 - m)))

Rule 2657

Int[(cos[(e_.) + (f_.)*(x_)]*(b_.))^(n_)*((a_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Simp[b^(2*IntPart[
(n - 1)/2] + 1)*(b*Cos[e + f*x])^(2*FracPart[(n - 1)/2])*((a*Sin[e + f*x])^(m + 1)/(a*f*(m + 1)*(Cos[e + f*x]^
2)^FracPart[(n - 1)/2]))*Hypergeometric2F1[(1 + m)/2, (1 - n)/2, (3 + m)/2, Sin[e + f*x]^2], x] /; FreeQ[{a, b
, e, f, m, n}, x]

Rule 4395

Int[((f_.)*sin[(a_.) + (b_.)*(x_)])^(n_.)*((g_.)*sin[(c_.) + (d_.)*(x_)])^(p_), x_Symbol] :> Dist[(g*Sin[c + d
*x])^p/(Cos[a + b*x]^p*(f*Sin[a + b*x])^p), Int[Cos[a + b*x]^p*(f*Sin[a + b*x])^(n + p), x], x] /; FreeQ[{a, b
, c, d, f, g, n, p}, x] && EqQ[b*c - a*d, 0] && EqQ[d/b, 2] &&  !IntegerQ[p]

Rubi steps \begin{align*} \text {integral}& = \left (\cos ^{-m}(a+b x) \sin ^{-m}(a+b x) \sin ^m(2 a+2 b x)\right ) \int \cos ^m(a+b x) \sin ^{-3+m}(a+b x) \, dx \\ & = -\frac {\cos ^2(a+b x)^{\frac {1-m}{2}} \csc ^2(a+b x) \operatorname {Hypergeometric2F1}\left (\frac {1-m}{2},\frac {1}{2} (-2+m),\frac {m}{2},\sin ^2(a+b x)\right ) \sec (a+b x) \sin ^m(2 a+2 b x)}{b (2-m)} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 6 vs. order 5 in optimal.

Time = 17.64 (sec) , antiderivative size = 2308, normalized size of antiderivative = 27.15 \[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\text {Result too large to show} \]

[In]

Integrate[Csc[a + b*x]^3*Sin[2*a + 2*b*x]^m,x]

[Out]

(AppellF1[-1 + m/2, -m, 2*m, m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Cos[(a + b*x)/2]^2*Cot[(a + b*x)/2]
^2*Sin[2*(a + b*x)]^m)/(2*b*(-2 + m)*(2*(AppellF1[m/2, 1 - m, 2*m, 1 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)
/2]^2] + 2*AppellF1[m/2, -m, 1 + 2*m, 1 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2])*(-1 + Cos[a + b*x]) +
 AppellF1[-1 + m/2, -m, 2*m, m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*(1 + Cos[a + b*x]))) + ((4 + m)*App
ellF1[1 + m/2, -m, 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sec[a + b*x]*Sin[(a + b*x)/2]^2*Sin[
2*(a + b*x)]^m)/(2*b*(2 + m)*((4 + m)*AppellF1[1 + m/2, -m, 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2
]^2]*(1 + Sec[a + b*x]) - 4*m*(AppellF1[2 + m/2, 1 - m, 2*m, 3 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]
 + 2*AppellF1[2 + m/2, -m, 1 + 2*m, 3 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2])*Sec[a + b*x]*Sin[(a + b
*x)/2]^2)) + ((4 + m)*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sin[a +
 b*x]^2*Sin[2*(a + b*x)]^m)/(4*b*(2 + m)*(2*(m*AppellF1[2 + m/2, 1 - m, 1 + 2*m, 3 + m/2, Tan[(a + b*x)/2]^2,
-Tan[(a + b*x)/2]^2] + (1 + 2*m)*AppellF1[2 + m/2, -m, 2 + 2*m, 3 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]
^2])*(-1 + Cos[a + b*x]) + (4 + m)*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/
2]^2]*(1 + Cos[a + b*x]))) + (2^(-3 + m)*Cot[(a + b*x)/2]*(Sec[(a + b*x)/2]^2)^(2*m)*(Cos[(a + b*x)/2]*(-Sin[(
a + b*x)/2] + Sin[(3*(a + b*x))/2]))^m*Sin[2*(a + b*x)]^m*((2 + m)*AppellF1[m/2, -m, 2*m, 1 + m/2, Tan[(a + b*
x)/2]^2, -Tan[(a + b*x)/2]^2] - m*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2
]^2]*Tan[(a + b*x)/2]^2))/(b*m*(2 + m)*(Cos[a + b*x]*Sec[(a + b*x)/2]^2)^m*((2^(-1 + m)*(Sec[(a + b*x)/2]^2)^(
2*m)*(Cos[(a + b*x)/2]*(-Sin[(a + b*x)/2] + Sin[(3*(a + b*x))/2]))^(-1 + m)*(Cos[(a + b*x)/2]*(-1/2*Cos[(a + b
*x)/2] + (3*Cos[(3*(a + b*x))/2])/2) - (Sin[(a + b*x)/2]*(-Sin[(a + b*x)/2] + Sin[(3*(a + b*x))/2]))/2)*((2 +
m)*AppellF1[m/2, -m, 2*m, 1 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2] - m*AppellF1[1 + m/2, -m, 1 + 2*m,
 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Tan[(a + b*x)/2]^2))/((2 + m)*(Cos[a + b*x]*Sec[(a + b*x)/2
]^2)^m) + (2^m*(Sec[(a + b*x)/2]^2)^(2*m)*(Cos[(a + b*x)/2]*(-Sin[(a + b*x)/2] + Sin[(3*(a + b*x))/2]))^m*Tan[
(a + b*x)/2]*((2 + m)*AppellF1[m/2, -m, 2*m, 1 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2] - m*AppellF1[1
+ m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Tan[(a + b*x)/2]^2))/((2 + m)*(Cos[a + b
*x]*Sec[(a + b*x)/2]^2)^m) - (2^(-1 + m)*(Sec[(a + b*x)/2]^2)^(2*m)*(Cos[a + b*x]*Sec[(a + b*x)/2]^2)^(-1 - m)
*(Cos[(a + b*x)/2]*(-Sin[(a + b*x)/2] + Sin[(3*(a + b*x))/2]))^m*(-(Sec[(a + b*x)/2]^2*Sin[a + b*x]) + Cos[a +
 b*x]*Sec[(a + b*x)/2]^2*Tan[(a + b*x)/2])*((2 + m)*AppellF1[m/2, -m, 2*m, 1 + m/2, Tan[(a + b*x)/2]^2, -Tan[(
a + b*x)/2]^2] - m*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Tan[(a + b
*x)/2]^2))/(2 + m) + (2^(-1 + m)*(Sec[(a + b*x)/2]^2)^(2*m)*(Cos[(a + b*x)/2]*(-Sin[(a + b*x)/2] + Sin[(3*(a +
 b*x))/2]))^m*(-(m*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sec[(a + b
*x)/2]^2*Tan[(a + b*x)/2]) + (2 + m)*(-1/2*(m^2*AppellF1[1 + m/2, 1 - m, 2*m, 2 + m/2, Tan[(a + b*x)/2]^2, -Ta
n[(a + b*x)/2]^2]*Sec[(a + b*x)/2]^2*Tan[(a + b*x)/2])/(1 + m/2) - (m^2*AppellF1[1 + m/2, -m, 1 + 2*m, 2 + m/2
, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sec[(a + b*x)/2]^2*Tan[(a + b*x)/2])/(1 + m/2)) - m*Tan[(a + b*x)/2
]^2*(-(((1 + m/2)*m*AppellF1[2 + m/2, 1 - m, 1 + 2*m, 3 + m/2, Tan[(a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sec[(a
 + b*x)/2]^2*Tan[(a + b*x)/2])/(2 + m/2)) - ((1 + m/2)*(1 + 2*m)*AppellF1[2 + m/2, -m, 2 + 2*m, 3 + m/2, Tan[(
a + b*x)/2]^2, -Tan[(a + b*x)/2]^2]*Sec[(a + b*x)/2]^2*Tan[(a + b*x)/2])/(2 + m/2))))/(m*(2 + m)*(Cos[a + b*x]
*Sec[(a + b*x)/2]^2)^m)))

Maple [F]

\[\int \csc \left (x b +a \right )^{3} \sin \left (2 x b +2 a \right )^{m}d x\]

[In]

int(csc(b*x+a)^3*sin(2*b*x+2*a)^m,x)

[Out]

int(csc(b*x+a)^3*sin(2*b*x+2*a)^m,x)

Fricas [F]

\[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\int { \sin \left (2 \, b x + 2 \, a\right )^{m} \csc \left (b x + a\right )^{3} \,d x } \]

[In]

integrate(csc(b*x+a)^3*sin(2*b*x+2*a)^m,x, algorithm="fricas")

[Out]

integral(sin(2*b*x + 2*a)^m*csc(b*x + a)^3, x)

Sympy [F]

\[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\int \sin ^{m}{\left (2 a + 2 b x \right )} \csc ^{3}{\left (a + b x \right )}\, dx \]

[In]

integrate(csc(b*x+a)**3*sin(2*b*x+2*a)**m,x)

[Out]

Integral(sin(2*a + 2*b*x)**m*csc(a + b*x)**3, x)

Maxima [F]

\[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\int { \sin \left (2 \, b x + 2 \, a\right )^{m} \csc \left (b x + a\right )^{3} \,d x } \]

[In]

integrate(csc(b*x+a)^3*sin(2*b*x+2*a)^m,x, algorithm="maxima")

[Out]

integrate(sin(2*b*x + 2*a)^m*csc(b*x + a)^3, x)

Giac [F]

\[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\int { \sin \left (2 \, b x + 2 \, a\right )^{m} \csc \left (b x + a\right )^{3} \,d x } \]

[In]

integrate(csc(b*x+a)^3*sin(2*b*x+2*a)^m,x, algorithm="giac")

[Out]

integrate(sin(2*b*x + 2*a)^m*csc(b*x + a)^3, x)

Mupad [F(-1)]

Timed out. \[ \int \csc ^3(a+b x) \sin ^m(2 a+2 b x) \, dx=\int \frac {{\sin \left (2\,a+2\,b\,x\right )}^m}{{\sin \left (a+b\,x\right )}^3} \,d x \]

[In]

int(sin(2*a + 2*b*x)^m/sin(a + b*x)^3,x)

[Out]

int(sin(2*a + 2*b*x)^m/sin(a + b*x)^3, x)